Explicit Euler method for solving the time dependent Schrodinger equation 
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Using an explicit Euler substitution, a system of differential equations was obtained, which can be 
used to find the solution of the time-dependent f-dimensional Schrodinger equation for a general 
form of the time-dependent potential. 
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In the framework of fundamentals of Quantum me- 
chanics, the unitary evolution given by the Schrodinger 
equation is conceivable only for quasi-isolated quantum 
systems, i.e. systems interacting via "classical" fields, 
which may be obtained using quantification procedures. 
Actually, these fields are time-independent, or have a 
harmonic time-dependence. The evolution of a quan- 
tum system interacting with a non-stationary environ- 
ment is given by master-type equations rather than time- 
dependent Schrodinger equations. However, the latter 
are often used in many practical applications from molec- 
ular physics, quantum chemistry, quantum optics, solid 
state physics, etc. 

In the recent paper [Q, I. Guedes obtained the ex- 
act Schrodinger wave function for a particle in a time- 
dependent 1-dimensional linear potential energy. For a 
Hamiltonian: 
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if Schrodinger equation is tested with the trial function: 

= Are[''(*)-^+''(*)l (2) 

one obtains two differential equations for the functions 
ri{t) and so one can find the solution after 
solving these equations. To this end, one has to find the 
initial conditions, say from ^(a:,0). However, \['(a;,0) is 
not, generally, of the type (||), so one has to decompose 
^(x, 0) using Fourier series: 
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next find the initial conditions corresponding to the gen- 
eral component tp{k) ■ e~* '^'^, then one has to solve the 
differential equation for ri{t,k) and ii{t,k), and finally 
compute the inverse Fourier Transform. 

Here I present an algorithm which can be used in solv- 
ing the Schrodinger equation for a general form of the 
time and position dependence of potential energy, with- 
out referring to the Fourier Transform. One can test 
Schrodinger equation using the following explicit Euler 
substitution: 


^'(x, t) = N ■ exp[Y^ a„(t) • a;"] 


I used the series expansion of the potential energy w.r.t. 
x: 


The initial conditions are given by: 

oo 

^(x, 0) = N ■ exp\^ a„(0) • x"] 
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Expanding in power series, Schrodinger equation gives: 
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[(n + 2) • (n + 1) • a„+2 + (7) 


+ V(fc + l)-(n-fc + l)-afc+i-a„_fc+i] + -V;(")(0,t) = 

k=0 

One can easily note that if one has non-zero initial coef- 
ficients only for n < uq (which includes :^Vx'"^\o,t)), all 
coefficients for n < 2 -no have to remain zero at every sub- 
sequent moment (one can perform any-order time deriva- 
tions in (0) at i = and obtain null values). Relation (0) 
can be very useful, either for analytical calculations, or 
for numerical (finite difference) calculus algorithms: 

a«,p+i = an,p + ^-7r—-^[{n + 2) • (n -f- 1) ■ a„+2,p + (8) 
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+ Y,(k+l)-{n-k+l)-ak+i.p-an-k+i.p] + —,Vj"'>{0,p-tQ) 

k=Q 

In (H) a finite time step to was chosen . It is obvious 
that one can iterate (^ in order to estimate {a„^p}„ for 
every moment p - to as functions of the initial coefficients 
{an,o}n given by (||). 


(4) [1] I. Guedes, Phys. Rev. A, 63, 034102 (2001). 


n=0 


1 


